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Abstract 

We investigate the low energy phenomenology of the lighter pseudoscalar ^4 5 in 
the NMSSM. The A\ mass can naturally be small due to a global U(1)r symmetry 
of the Higgs potential, which is only broken by trilinear soft terms. The A® mass 
is further protected from renormalization group effects in the large tan f3 limit. 
We calculate the b — > sA® amplitude at leading order in tan/3 and work out the 
contributions to rare K, B and radiative T-decays and B — B mixing. We obtain 
constraints on the A® mass and couplings and show that masses down to O(10) 
MeV are allowed. The 6-physics phenomenology of the NMSSM differs from the 
MSSM in the appearance of sizeable renormalization effects from neutral Higgses to 
the photon and gluon dipole operators and the breakdown of the MSSM correlation 
between the B s — ► fJ. + n~ branching ratio and B s — B s mixing. For A\ masses 
above the tau threshold the A® can be searched for in b — > st + t~ processes with 
branching ratios < 10 -3 . 
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1 Introduction 



Sizeable flavor changing neutral current (FCNC) effects in meson decays arise in the 
minimal supersymmetric Standard Model (MSSM) at large tan/3, e.g. PQ-[E]- In this 
model, the amplitude of exchanging neutral Higgses between down-type fermions /, 
i.e. down-type quarks or charged leptons 

J2 (9ffs) 2 oc cos 2 (/3 - a) sm 2 (/3 - a) | 1 _ Q ^ 

S=h°,H°,A° m S m h° m H° m A" 

vanishes. Here, m s , gjf S denote the Higgs masses and couplings to a fermion pair, 
respectively and a is the scalar mixing angle. Eq. ([TJ implies that the Wilson coefficients 
for b — > s£ + £~ decays from scalar and pseudoscalar boson exchange in the MSSM at large 
tan/3 are equal with opposite sign 0,0 • If the relation is broken, interesting effects via 
operator mixing are induced In particular, the dipole operators responsible for 
b — > S7 and b — > sg decays receive sizeable contributions from the neutral Higgs bosons. 
Furthermore, specific contributions to B — B mixing from scalar exchange arise. This 
happens in the presence of more Higgses, such as in the next-to-minimal supersymmetric 
Standard Model (NMSSM). 

The NMSSM is the MSSM extended by a singlet N, with the superpotential [HJ E! 

W = QY U H U U + QY d H d D + LY e H d E + \H d H u N - ^kN 3 (2) 

The physical NMSSM Higgs sector consists of three scalars h°, H® 2 and two pseudoscalars 
A\ 2 - As in the minimal model, tan/3 = v u /v d denotes the ratio of Higgs doublet vevs 
v u =< H® >= v sin/3 and v d =< H® >= vcos/3, where v = \/2mw / g — 174 GeV. The 
Higgs potential 

Vhiggs = V so ft + Vf + Vd (3) 

where 

Vsoft = m 2 Hd \H d \ 2 + m 2 Hu \H u \ 2 + m 2 N \N\ 2 - (XA x H d H u N + h.c.)- (hA k N 3 + h.c.) (4) 
V F = \X\ 2 (\H d \ 2 + \H u \ 2 ) \N\ 2 + \XH d H u - kN 2 \ 2 (5) 

Vo = (\H d \ 2 - \H U \ 2 ) + 9 ^\HlH d \ 2 (6) 

has a global U(1)r symmetry in the limit of vanishing soft terms A^, A\ — > ^U]. If this 
symmetry is broken only slightly, the model naturally contains a light pseudoscalar. Its 
mass is given as 

m 2 A0 = 3kxA k + 0{——) (7) 
1 tan p 

where x =< N > denotes the vev of the singlet. Note that a small A k remains small 
under renormalization group running and thus protects m^o. 
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Lower bounds on CP-odd scalar masses are not very stringent and can be as low 
as ~ 100 MeV jllj . Since the coupling h?A\A\ is not suppressed the scalar Higgs 
predominantly decays into the lighter pseudoscalars. This has important consequences 
for the Tevatron and LHC Higgs searches [TTH I12j. 

The motivation for this work is to find out how and to what extend the NMSSM would 
signal itself in rare 6-decays and at the same time, whether existing data provide already 
bounds on the NMSSM parameter space. We employ the large tan/5 > 30 and small 
Ak <C my/ , x limit and no flavor or CP violation other than in the CKM matrix ( "minimal 
flavor violation"). Since a small A\ is not stable under radiative corrections, we do not 
expand in small A\ and keep it finite. Our study is based on mostly generic features of 
the NMSSM. Specific analyses of the NMSSM particle spectrum and parameter space 
have been carried out in a GUT framework [T3] at large tan (3 [Hj , with gauge mediated 
SUSY breaking and with anomaly mediation [T^]. For Higgs production in rare 
6-decays in other models, see e.g. [T7| ITS]. 

This paper is organized as follows: In Section |2] we calculate the amplitude for 
b — ► sA\ decays at large tan f3. We discuss the NMSSM parameter space in Section El 
Phenomenological bounds from FCNC decays, B — B mixing and T-decays are worked 
out in Section |3J In Section |5] we investigate the impact on semileptonic and radiative 
rare 6-decays. We also analyse how much the MSSM tree level relation Eq. is broken 
by loop corrections. We conclude in Sectional Feynman rules and the NMSSM particle 
spectrum at large tan /3 and auxiliary functions are given in Appendix El and El In 
Appendix we give decay rates of the A® and 6-decay branching ratios. 



2 The b —> sAi amplitude at large tan (3 

The amplitude for a FCNC b — > s transition into the lightest CP-odd scalar A\ in the 
NMSSM is induced at one-loop. In the large tan /3 limit, only two diagrams remain to be 
calculated, which are shown in Figure ^ (We neglect the strange quark mass). Feynman 
rules are given in Appendix |A~TJ see also HHj for the MSSM and [H] for the NMSSM. 

The stop chargino wave function correction is identical to the corresponding one in 
the MSSM. Since the coupling of the A\ to down-type fermions is order (tan /3)°, the 1PR 
diagram contributes to the b —>■ sA® amplitude at order tan/3. The vertex correction 
shown in Figure^is the only 1PI diagram linear in tan j3 because % the H ± W^A\ coupling 
is 1/ tan/3 suppressed since the A\ is predominantly the gauge singlet (the H + H~A®, 
W + W~A® vertices are forbidden by CP), ii the coupling of the A\ to up-type quarks 
is 1/ tan 2 /3, Hi the coupling of the A\ to up-type squarks is 1/ tan/? which can be seen 
from the F-term contribution \dW/dH u \ 2 and iv the only tan j3 enhancement comes from 
the bRtiHd or bntiHa vertices. 

We obtain the following amplitude 

iA{b -> sAf) = -A%V tb V: s C A -^—s L b R A\ (8) 
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Figure 1: The leading b — > sA® diagrams at large tan/3 in the NMSSM. 
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.tan f3m,b 



— ^2 Xi + Xm t sin Of cos Of Y\j 

X »=1,2 i,j=l,2 



(9) 



and 5_ parametrizes the A\bb coupling, see Eq. (|A-4j) . The X, Y terms in Eq. (jHJ) result 
from the wave function and vertex correction, respectively . They are written as 

Xi = m Xt U i2 V2m w Vn ^-D 3 {y Cyi ) + £> 3 (y 1(i ) cos 2 { + D 3 (y 24 ) ^ 0f) 

- mtV^sin^cos^ p 3 (yii) - Ajfe))] (10) 
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and mf , mj, m X; denote the stop, scharm and chargino masses. The stop mixing angle Of, 
the chargino mixing matrices Z7, 1/ and the loop functions .D2 > -D3 are defined in Appendix 
IHl We used unitarity of the CKM matrix and neglected squark mixing other than for 
stops and mass splitting between the first two generations. The b — > dA\ amplitude is 
obtained by replacing V by l d' everywhere in Eq. (JHJ). The s — > cM^ amplitude is given 
correspondingly with also changing m& to m s in Eq. (0). Note that our calculation holds 
for \6±v/x\ < tan/?, see Appendix lAl E.g. for larger values of 5_v/x the A\ looses its 
mostly-singlet nature and more b — > sA® diagrams need to be calculated. 

The coupling Ca vanishes if the super GIM mechanism is active, that is either all 
squark masses are degenerate or ma = and Of = (or it) or 
We estimate the generic size of Ca with order one stop mixing as 



m f ~ 2 and Of = n/2. 



\Ca\ — O (5_ tan f3m h m t — -) + 0(X tan (3m b m t ) 



x 



(13) 
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Since Ax is the NMSSM /i-term which sets the mass scale for the charginos, both terms 
are of comparable size. 



3 Viable points in the NMSSM parameter space 

The relevant NMSSM parameter space consists of A, k from the superpotential, the soft 
breaking terms A\, Ak, the gaugino mass m 2 , stop and scharm masses, the stop mixing 
angle Qi and tan (3. We evaluate all parameters at the electroweak scale. 

The dimensionless couplings A and k run towards smaller values, e.g. A 2 + k 2 < 0.6 
at the electroweak scale for A, k < 2tt at the high, GUT scale [21]. We use |A|, \k\ < 1. 
Similar to the MSSM, electroweak symmetry breaking at large tan j3 requires 

m 2 Hu = -(Ax) 2 - ^ (14) 

and therefore the product of A and x should not exceed ~ 0(1) TeV to avoid fine tuning. 
On the other hand, the chargino mass scale is driven by Ax, which should be at least 
(9(100) GeV by experimental search limits. We assume the singlet vev x to be of the 
order of the Fermi scale v , or at least not smaller than 100 GeV and not bigger than 3 
TeV. If x exceeds this value, its relation to the other vevs becomes unnatural and the 
model does not give a solution to the fi problem [H]. Hence, the size of A is bounded 
from below as |A| ^ few • 10 -2 |TH| I21j. Further, the extremization condition 

< = -AV + ^) + ^ + f (15) 

where we defined 

m 2 A = \(A\ + /cx)xtan/3 (16) 
implies some cancellation among the tan (3 enhanced terms as [T4] 

0(100 - 1000 GeV) 
A x + kx i 17 

tanp 

Note that m A sets the scale for the heavy Higgses A\, H® and H ± , see Appendix 1X1 

The NMSSM is further constrained by non-observation of Higgses and superpartners. 
At large tan (3, the mass of the lightest scalar at tree level is given as 

2 2 /.on 

m h0 ~m z - — (18) 

where we expanded Eq. ()A-7|) in m|/4A; 2 x 2 1 and X 2 v/k 2 x -C 1. In this approximation 
also m 2 H0 ~ 4A; 2 x 2 + and the scalar mixing ang le 9 is small. Like in the MSSM, the h° 
tree level mass cannot be bigger than the Z-mass because the raising of its upper bound 
in the NMSSM is suppressed by large tan (3. To be phenomenologically viable, myfi has 
to be lifted by radiative corrections above the current search limit as in the MSSM [22\. 
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We require the scalar tree level mass to be bigger than 89 GeV, which favors small A or 
X/k less than one. We allow for |m 2 | < 1 TeV and check that the charginos are heavier 
than 90 GeV. We treat the pseudoscalar masses m^o and m A o with m A o > 130 GeV as 
free parameters, i.e. adjust Ak and A\ accordingly. The squark masses and stop mixing 
angle are effective parameters with > 90 GeV and m t ~ 2 , mg ~ 1 TeV and we do not 
relate them to fundamental parameters in the Lagrangian. 

The down-type fermion-A^ vertex is proportional to 8_v/x, see Eq. flA-9|) . From 
Eqs. (ITB|) and (|T7jl we obtain 

v v 

X X 

where the second equation is a good approximation for not too large m A < 500 GeV. It 
then gives a lower bound on \S-.v/x\. In particular, for tan/3 = 30, m A ~ 500 (200) (130) 
GeV and the ranges of parameters given in the preceding paragraphs, we obtain \5~v/x\ > 

0. 1 (1) (3). For larger values of m A cancellations between the two terms in Eq. (JT§|) are 
possible. Note that the tan/3 factor is only a formal enhancement, since it is cancelled 
by the one in m\. We find that \S^v/x\ < 62 (16) for m A > 500 (1000) GeV. Note that 
the small A\ <C kx limit with 5_ ~ —2 makes its hard to satisfy Eq. (jl7j) . 

4 Phenomenology of the light A\ 

We work out constraints on the mass of the A\ in the NMSSM at large tan/5 from A\ 
production in rare decays ( Section l4.1j) . B — B mixing ( Section l4.2j) and B s — > fi + 
decays fSection l4.3j) . We make use of the b — > sA® amplitude calculated in Section |21 We 
scan the parameter space in the regions discussed in Section El All FCNC bounds can 
be evaded by a sufficiently tuned-in super GIM mechanism, see Section To quantify 
this, we demand in our numerical analysis for the mass splitting m~ t2 — > 50 GeV 
while varying and for the stop mixing e<0j<7r/2 — eor7r/2 + e<%<7r — e with 
e = 0.05. Bounds from other processes are discussed in Section l4~4l 

Many experimental constraints we use here apply only if the A\ is sufficiently stable, 

1. e. leaves the detector as missing energy. This happens if the pseudoscalar width is 
smaller than E A / (m A o d), where d ~ O(10) m is the size of the detector and E A the A\ 
energy in the lab frame. We work out bounds on m A o as a function of \5-.v/x\. If this 
coupling gets smaller, the pseudoscalar decay rate decreases, and a heavier Higgs will 
become missing energy and vice versa. For decay rates of the A®, see Appendix O 

For Higgs masses below 2m M only the e + e~ and 77 decay channels are relevant. 
(The A\ — > 7r°7 decay is forbidden by CP and angular momentum conservation and the 
A\ — > 7r°77 decay is suppressed with respect to the dielectron mode by phase space and 
powers of a). The 77 mode can compete with A\ — > // decays only near the dimuon 
threshold. This weakens the missing energy bounds in that region. 

The point we want to make is to show that in the NMSSM A\ masses in the GeV 
range and below are not ruled out. This is summarized in Figure 121 For details see 



k\x 2 
' — T 



tan p + 1 



±3 ^tan/3 

m A 



(19) 
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the following subsections. All experminental bounds are taken at 90 % confidence level. 
The requisite b — > A\ branching ratios are given in Appendix O We recall that our 
approximation breaks down if \8±v/x\ approaches tan/3. 
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Figure 2: Constraints on the A\ mass as a function of \8_v/x\ at tan/3 = 30 in the 
NMSSM. Shaded regions are excluded. The left bottom corner is excluded by rare K- 
decay, see Eq. \21\) . The triangular region to the lower right is obtained from radiative 
T(ls) decays, see Eq. \35)) . The region to the left of the vertical dashed blue lines can 
only be reached if is bigger than the value indicated, see Section^ Constraints from 
Arrid are given for > 500 GeV and tua > 1000 GeV. We also show the missing 
energy condition for B —> K decays given in Eq. (H?J) (dashed green line). The vertical 
dashed lines indicate m A o = Im^ and 3m n . 

4.1 Rare K and 5-decays 

If the Higgs boson is light enough, it can be produced in b — > sA® or s — > dA\ processes. 
We analyze what bounds exist depending on the mass of the A\. 
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4.1.1 2m e < m A o < 2m ll 

When produced in rare B-meson decays the A\ decays outside of the detector if 

m A o < 17 MeV/\5_v/x\ (20) 

In this region the CLEO bound B(B -> K%X°) < 5.3 • 10~ 5 [21] applies. There is a 
similar missing energy bound from BaBar B(B~ — > K~vv) < 7.0 • 10~ 5 1 . We find 
that masses in the range given in Eq. (J20)) are disfavored since the B — > KA\ decay, see 
Eq. (jC-5|) . would happen too rapidly for most of the parameter space, although cannot 
rigorously be excluded. We stress that the size of the coupling C A can be quite large, 
see Eq. (JI3J) and already B(B -> X S A\) < 1 cuts out a fraction of NMSSM points. 

Rare decays into e + e~ constrain Higgs masses below the muon threshold. However, 
the measurements of the inclusive B — > X s e + e~ branching ratios contain cuts on the 
dilepton mass m ee > 2m ii [23123- I n the analysis of B — > K^e + e~ decays Belle applies 
m ee > 0.14 GeV [2Zj, whereas BaBar j2H] has no cut, but the efficiency is low in that 
region due to conversion photons. Likewise, measurements of K + — > 7r + e + e~ decays 
employ a high mass trigger |2H|. Since also close to 2m M the two-photon decay of the A\ 
becomes sizeable, we do not take the e + e~ data into account. 

The bound B(K + -> tt + A\) < 4.5 ■ 10" 11 [SU] is applicable if the A\ becomes suffi- 
ciently stable to escape the detector. This happens for masses 

m A o Z 5MeV/\5^v/x\ (21) 

which then are excluded. The i^-decay bound is five orders of magnitude better than the 
one from B — > K decays, because the CKM and mass suppression of the K — > tiA\ decay 
rate is compensated by the difference in life time \V t d/V ts \ 2 (mK /^_b) 3 t(-^ + ) l T {B + ) — 
0.24 :3I], see Eq. ()C-5|) and its K — ► 7r counterpart. 

4.1.2 2m M < m A Q < 2m T 

A\ decays into a muon pair are included in B — > X s fi + fi~ signals. Comparison of the 
B -> X S A\ branching ratio, see Eq. (jC-ijl . with the data B(B -> X e fi + fi~) < 10.4 • 10" 6 
[3 123 I2H] shows that this is very unlikely. The same happens in K —>■ decays, 
which for m A o < — can hide a pseudoscalar decaying into muons. With B{K + — ► 
7t + {i + fi~) < 10.4 • 10~ 8 [31] only a tiny number of points survives the scan. All allowed 
points are at the GIM boundary % ~ 7r/2, which is set by our value of the cut-off e. 

Above the 3ti threshold sizeable hadronic decays open up. (The A\ — > 2n^ decay is 
suppressed with respect to the dimuon channel by phase space and a, whereas A\ —>■ 2n 
decay is forbidden by CP invariance.) For the A\ decaying hadronically into a strange 
final state we use B(b — ► sg) < 9% [S2]- This thins out the NMSSM model space for 
3m n < m A o < 2m T , but cannot exclude this region. (We use B{b sA?) > Bib — > dAf).) 

lr rhe experimental cut on the K momentum \pk\ > 1.5 GeV is no restriction for light m A a <c tub 
discussed here. 
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4.1.3 2m T < m A o < m# 

If the A\ is above the tau threshold, most of the time it decays into t + t~ because its 
coupling to cc is tan 2 ft suppressed. Similar to the constraint on the hadronically decaying 
pseudoscalar, see Section BjH2| the mildly model-dependent bound B{B — > X s t + t~) < 5 
% [33] is not a challenge to the light CP-odd Higgs scenario. 



4.2 NMSSM neutral Higgs contributions to B — B mixing 

We calculate the contribution to B — B mixing from pseudoscalar A\ 2 and scalar h°, Hf 2 
Higgs exchange in the NMSSM at large tan ft. It arises at two-loop from double insertion 
of the FCNC s6-Higgs vertices such as generated by the diagrams in Figure for the 
A\ and an intermediate boson propagator. The dominant diagrams induced by the 
heavy Higgses, i.e. the ones other than the lightest CP-odd scalar are the wave function 
corrections contributions with A®, exchange, see the Feynman rules in Appendix lA.il 
They can compete with one-loop contributions such as the Standard Model (SM) box 
diagrams due to their tan 4 /3 enhancement. Contributions from h°, if° are subleading in 
tan ft. We use an effective Hamiltonian (q = d, s) 

^B^^^.^^ (22) 

i 

where some of the relevant operators are written as, see e.g. |H] 

Q VLL = (qLl,b L )(q L Yb L ) (23) 
Ql RR = {q L b R ){q L b R ) (24) 
Qf LR = (qRb L )(q L b R ) (25) 

The SM contribution is in the coefficient C VLL . The A®, H® masses are degenerate at 
large tan/3 and their respective contributions to Qf cancel each other just like in the 
MSSM, see Eq. (JTJ). They do, however, contribute to the operator Qf LR at order m q /rrib, 
and are important for S s -mesons. (This is the famous double penguin (DP) contribution 
of the MSSM 0E].) 

We obtain at order tan 2 ft/w? A o in the NMSSM from A\ boson exchange 

C! RR {^) = - j^-rV ( 26 ) 

An 2 m^m A o 

at the high, electroweak (matching) scale fi t . Finite widths effects are neglected. We 
define the size of the B q — B q mass difference Am, with respect to its SM value as 

{Am q )sM 
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where 

= & M <*> (28) 

and S (>t) = 2.38 and Pf LL = -0.37 [6 . In Eq. the NMSSM contribution to Am d 
by neutral Higgs exchange in the m q = limit has been given. To be in agreement 
with data we require fd > —0.6 (fd is negative). This includes 20 % uncertainty and 
allows for cancellations between the A\ contribution and the charged Higgs, chargino 
boxes and the double penguins. We assume similar sizes as in the MSSM, where —0.2 < 
fd* + fd + f d P ~ 0.4 jH]. We find constraints for larger values of \5-v/x\ and itla > 500 
GeV, which are displayed in Figure El for tan/? = 30. The other branch with 1 + fd < 0, 
where the NMSSM correction is larger than the SM box gives very similar constraints 
and is not shown. The leading A\ contribution to B — B mixing is universal in minimal 
flavor violation, f d = f s , since we neglect light quark masses. 



4.3 B s — ► fi + fi decays 

We work out the contributions to B s — > fi + n~ decays from neutral Higgs exchanges in 
the large tan (3 limit of the NMSSM. With the effective Hamiltonian 

% = -|^E^< ( 29 ) 



where 



e 2 _ , - nn „ e 2 



°s = T7T^^b R ££ , P = —s L b R lld (30) 
lbir z 167r / 

we obtain at the electroweak scale (in parentheses is given the particle that induces a 

particular Wilson coefficient) 



where 



C s = C S (H° 2 ) = -C P (A° 2 ) , C P = C P {A\) + C P (A° 2 ) (31) 



v ( vS 

Cp(A° 1 ) = rribmitan f3- 



x 



Am^ sin 2 8 W \ x ) m 2 Bs — m 2 A0 



^2 Xi + ^ m t sin 9i cos Q% ^ 

x i=l,2 z,i=i,2 



(32) 



C P {Al) = m bmi tan 3 B ^ £ X, (33) 

4m A m w sm v w i=12 

The expressions for X and K are given in Sectional Our result for the A 2 , H 2 contribu- 
tions agrees with the corresponding MSSM calculations Note that the contributions 
from A 2 and H 2 are equal with opposite sign. Similar to B — B mixing discussed in 
Section 14. 2\ the scalars h° and Hf contribute at subleading order in tan (3. 
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The coefficients Cs,p are model-independently constrained by data on the B s — > fi + {i 
branching ratio. With B(B S — > /i + /i~) < 5.8 • 1CT 7 [34J we obtain at the scale /i = mw 



C s \ 2 + \Cp + 5 w \ 2 < 1.3 



B{B S 


-> fi + n ) 


1/2 


"238 MeV" 


5.8 


x 10~ 7 







(34) 



Here, 5 W stems from the operator C 10 oc s if) JolQ-I* 1 see [7] for details. We find with 
<5f M = -0.095 2 at tan/3 = 30 the upper limits \6-v/x\ < 42 (16) for m A = 500 (1000) 
GeV, which are are weaker than the corresponding Arrid ones. The expressions for the 
CKM suppressed Bd — > decay are readily obtained. Its experimental constraint is 

not as good as the B s one, but we can cut out both m A o ~ m Bd and m A o ~mg s . 



4.4 Non-FCNC bounds 

The bounds from radiative T-decays apply if the A\ leaves the detector unseen |3*T| I3fij . 
Due to the larger boost the critical width to do so is larger than in 5-meson decays by 
m r /m B . We use B(T(ls) -> A° 7 ) < 1.3 • 10" 5 jSi and obtain with Eq. (|(T7jl 

[<y_w/aj| < 3.7 for m A o < 23 MeV/\S^v/x\ < 2m M (35) 

Furthermore, we get an upper bound \8-v/x\ ^ 100 from B(T(ls) — > A°7) < 1. 

Mass bounds from hadronic collisions are not better than few to 200 MeV and astro 
physics gives m A o > 0.2 MeV |31j . which contain some model dependence. 



5 Implications for b — > s£ + £ and b — > S7, ^ 

Similar to the operators Os,p discussed in Section 14.31 the NMSSM Higgs sector also 
induces contributions to 4- Fermi operators with quarks and leptons (/ denotes a fermion) 

0{ = s L b R f R f L , O f R = s L b R f L f R (36) 

where 

= ™ (37) 

These couplings arise in the NMSSM at large tan /?, where 

C s - C P = -2C P (A° 2 ) - C P (A1) , C s + C P = C P {A\) (38) 

This is different from the MSSM, where the A\ contribution is absent and the sum of 
Cs and Cp and hence C R vanish. We discuss corrections to this tree level statement in 
Section Um All Wilson coefficients refer to the Hamiltonian in Eq. (}2T?j) and are evaluated 
at the scale \i = mw unless otherwise stated. 

2 Susy effects are not tan/3 enhanced in Oiq and are small with minimal flavor violation |35j . 
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The constraint given in Eq. (JMj) implies for the Wilson coefficients for 6-quarks (we 
update the findings of Ref. [7] with the improved B s — > bound |34j) 

^\C h L ? + \C b R \ 2 < 0.03 (39) 

The operators O b L R enter radiative and semileptonic rare b — > s decays at one- loop [71l3*7j. 
With the bound in Eq. (|3*9*j) the new physics effect from O b L is small, at the percent level 
[7j. However, the renormalization effect induced at leading log by O b R can be large for 
the photon and gluon dipole operators Of and which can be written as 

Or = j^m.SLa^bnF^, O s = ^- 2 m b s La a, u T^b Rp G a ^ (40) 

To be specific, we normalize their coefficients to the ones in the SM, and denote this 
ratio by £, such that £ SM = 1. With (see jTj for details) 

£ 7 ( m& ) = 0.514 + 0.450 &{m w ) + 0.035 £ 8 (m w ) - 2.319 C b R , (41) 
£ 8 ( mfe ) = 0.542 + 0.458 &(m w ) + 19.790 C b R (42) 

and Eq. (J3T?|) corrections of up to 7% and 59% to £ 7 and £ 8 are possible. This has impact 
on the extraction of Wilson coefficients in b — > s^y, b — > sg and b — > s£ + £~ decays [7j. 
For a full analysis of these decays, also the matching contributions to CV^ from neutral 
Higgs loops in the large tan (3 NMSSM have to be calculated. Note that tan (3 enhanced 
corrections to the 6-quark mass, CKM elements and FCNCs from non-holomorphic terms 
arise [H El IE] . We leave this for future work. 



5.1 Estimates of Cs + C P and C h R 

We work out the NMSSM reach in Cs + Cp by taking into account all constraints 
discussed in the previous Sections 01 and |U The value of Cs + Cp can saturate its upper 
bound given in Eq. for large ranges of the parameter space. If the A\ gets very light, 
however, the b — > sA® coupling Ca has to decrease and Cs+Cp is small, e.g. for m A o = 10 
MeV is \Cs + Cp\ < 0.06. For intermediate masses the A\ contribution dominates over 
the one from the heavy pseudoscalar, that is \Cp\ ^> Cs and \C R \ ~ \C b L \ < 0.024. This is 
illustrated in Figure EH where we show Cs + C R as a function of \j\Cs\ 2 + \Cp + <5fo M | 2 oc 

JB(B s — > // + )tt _ ) , see Eq. (HHJ), for tua = 500 GeV and different values of m A o- 

In the MSSM the size of Cs + Cp is driven by the relation Eq. which is not pro- 
tected from radiative corrections. To study their size we employ the two-loop calculation 
encoded in FeynHiggs v. 2.02 [SHj. By scanning the MSSM parameter space we find 

l ^ 3 + ^ P \mssm < 0.2 or \C s + C p \mssm < 0.08 and \C b R \ MSS M < 1.3 • 10~ 3 (43) 

The smallness of Cs + Cp is a feature of the Higgs sector of the MSSM. It holds also with 
flavor violation beyond the CKM matrix. As a result, the logarithmic renormalization 
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Figure 3: The correlation between C s + C P and ^\C S \ 2 + \Cp + 5f^\ 2 in the NMSSM 
for tan/? = 30, = 500 GeV and m A o = 0.1,1 and 10 CeK Also shown is the 
experimental upper bound given in Eq. (dashed line). 

of the dipole operators from neutral (pseudo)scalars is tiny in this model. For example, 
consider additional right handed currents, which induce contributions to the helicity 
flipped operators 0' i: i.e. the ones obtained from Oj with right R and left L chiralities 
interchanged. In this case, C' s — C' P mixes onto the flipped dipole operators 0' 7 8 , but 
C' s = C'p in the large tan p MSSM |S] . 

6 Conclusions 

We investigated the phenomenology of the light pseudoscalar A\ which lives in the 
NMSSM spectrum at large tan (3. The A® has suppressed gauge interactions but couples 
to Higgses and down-type matter. We calculated the b — » sA® amplitude at leading 
order in tan f3. Based on this, we estimated the NMSSM contributions to rare K, B and 
radiative T-decays with the A\ in the final state, B s — > p^yT decays and B — B mixing. 
We showed that low energy data provide constraints on the A\ mass and couplings, but 
leave masses down to (9(10) MeV viable, see Figure EJ The A\ predominantly decays 
into r + r~ for 2m T < m A o < 2mb, light hadrons for 3m n < m A o < 2m T and e + e~or 77 for 
2m e < m A o < 2m M . In the latter range, the A\ can live long enough to leave detectors 
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undecayed, depending on S^v/x. For masses within 2m M < m A o < 3m,r the A\ decays 
mostly into muon pairs. Like the one from B — > K decays given in Eq. (j2*U)l . this mass 
range has very tight FCNC constraints, see Section Ell. 2| but is not ruled out. 

The A\ can be searched for with improved measurements of T-decays or B — > K plus 
missing energy. The latter needs a high A-momentum cut to suppress the B — > Kvv 
background. For m A o above the *f>' mass the pseudoscalar can be seen in b — > st + t~ 
decays. The required sensitivity for e.g. the B — > -A s t + t~ branching ratio is B(B 
X s n+n-)m 2 T lm\ < 10" 3 . 

The NMSSM has different implications for 6-physics than the MSSM. In particular, 
the leading log neutral Higgs contribution to radiative b — > S7 and b — > sg decays is 
tiny in the latter, but can reach experimental upper limits in the former, see Section 
15.11 Furthermore, the MSSM correlation between the B s — > fi + [i~ branching ratio and 
B s — B s mixing [SH] breaks down due to the additional pseudoscalar. For example, for 
small \Cs/Cp\ the lighter CP-odd Higgs dominates the B s — > rate, which can be 

anything up to the experimental bound, see Figure El At the same time Am s is near its 
SM value because the leading A\ contribution is independent of the light quark flavor and 
constrained by Am^, and the double penguin from A^ is suppressed. This is in contrast 
to the MSSM, where a SM-like Am s implies an upper bound on B(B S — > 
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port and Bogdan Dobrescu for collaboration at an early stage of this project. G.H. grate- 
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A Higgs spectrum and couplings 

We give the tree level Higgs spectrum, mixing angles in the minimal flavor and CP vio- 
lating NMSSM at large tan (3 and A k <C m w , x. The mass matrices in gauge eigenstates 
can be seen in j^]. 

The mass eigenstates of the pseudoscalar mixing matrix can be written as 



A\ \ I cos 7 sin 7 \ A 
Al) \- sin 7 cos7 / I Aj 



(A-1) 



where Aj = ImN/v^, ^4° = v^sin (3 Im H^ + cos/?ImH°) and the Goldstone boson is 
given as G° = v^eos (3 Im H[} — sin/JImHjJ). The mixing angle and masses read as 



a u 1 u 

7 = — + ^ - + 0{- — 2^5) i.e. shi7~l, cos 7 ~ — ^ Q 5- (A-2) 



71 + V 5 + 0( ^ 

2 xtan/3 tan 2 /?' xtan/3 

and 

m A o = 3kxA k , m A o = m A (A-3) 

where we defined 

^ A x + kx v ' 
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and m\ is given in Eq. (|16j). 

The scalar mass matrix can be diagonalized analytically in the large tan/5 limit by 
first decoupling the heaviest state and then rotating the remaining 2 by 2 block by the 
angle 9 along the lines of Ref . [25 • The result can be written as 



fh°\ 


/ 






U 2 v 


\ 



tan (3 



tan j3 



with the mixing angle and scalar masses 
4X 2 vx 



^(sin^ + -5+ cos 9) sin^ cos 9 



— 1 —V £ 

tan f3 x tan (3 + / 



V 



/ ReH° 

ReH° ln 
\ ReN-x / 



(A-5) 



tan 26> 



4Px 2 



TO; 



4JfeV + m| =F v / (4fc 2 a; 2 - m|) 2 + 16A 4 x 2 t> 



The mass of the charged Higgs is given as 

2 2,2 \ 2 2 

m#± = m A + m w - A v 

A.l Feynman rules 



V4^2„,2 



(A-6) 
(A-7) 

(A-8) 



Feynman rules can be read off the Lagrangians given at leading order in tan/3. Note 
that A® ~ — A° MSSM in this limit. 

Couplings to up (u) and down (d) type fermions 



.gm u 1 / i 



C 



(h°,H?)dd 



2mvi/tan / 9 \xtan/3 



gm d 



5-A(, A u 2 u~f 5 u 



(A-9) 
(A-10) 



, (cos G - -6+ sin 9)h°, 
2m w \ x 



(sin 9 + -5 + cos 9)H^ tan /3F 2 U ) dd 



£(h°,F° 

Couplings to charginos 



gm u 



cos9h°,sm9H°, 



m 



(fcO^--— , i; 



MM 



(A-ll) 
(A-12) 



+i^A° 1 xf [U l2 V j2 L - U j2 V i2 R] x + 



(A-13) 



where L, i? = (1 =f 7s)/2 are chiral projectors. 
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B Conventions, loop functions 

The chargino mass matrix is written as ([imssm = —Xx) 

M ( m 2 V2m wS in[3\ 
\ \ / 2m w cos j3 —Xx I 

It is diagonalized by the orthogonal matrices U, V (we do not include beyond CKM CP 
violation) 

f/M x± \/ T = diag(m Xl ,m X2 ) (B-2) 
The stop mixing matrix is given as 



;b- 3 ) 



ti \ I cos 9$ sin 9f \ I t L 
t 2 ] l^-sin^- cos9- t ] \l R 

Here, £1,2 are the mass and Il,r the gauge eigenstates. 
The loop functions are defined as 

D ^yi = {1 -tC- V ) + ix " y) - D ^ = ~\ < B - 4 > 

D 3 (x) = ^, C 3 (l) = -1 (B-5) 

1 — X 

C Decay rates 

The rate of the light NMSSM pseudoscalar into down-type fermions is given as 



m) 



r(A^ff) = -^-5„) m Af n^-^.r (C-l) 

where r = 1 for leptons and r = Nc for quarks. The decay rate into up-quarks is 
1 / tan 4 f3 suppressed. The rate into two photons reads as 

rK-77) = ^^olE^I 2 (c-2) 

where for i = e,fx,r and d, s, b loops Jj = rQ^KjF(Kj)<5_t>/x, = mJJm 2 A0 and Qi is the 
charge of the fermion. The function F(k) can be seen in JUI]. It assumes the limits 

for k < 1 

kF{k) = for k > 1 (C-3) 
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Higgsino loops contribute as 1% = \^Ui2Vi2\mw/m x K xF{ K x)- ^ follows from Eq. ()C-3|) 
that near m A o < 2m^ the 77 rate is dominated by the muon loop. Contributions from 
up- type quarks are suppressed by 1/ tan 4 /?. 

The decay rates for inclusive and exclusive b — > sA® FCNCs read as 

T(B - X S A X ) - 2l07r 5 \ C a\ —3 (C- 4 ) 



Y{B - KA\) = ^ 6 y J V ^l 2 ^l/o(ml)l 2 ( m| ~ m ^ ) 2 
2 w 7i° m B nib 



(C-5) 



where the form factor / parametrizes the matrix element 

< K(p K )\s L b R \B(p B ) >= l( !^^k )f ((p B - PK f) (C-6) 

2 nib 

Here, pk denotes the three momentum of the Kaon and /o(0) ~ 0.3 to 0.4 
The branching ratio for radiative T decays is given as, e.g. [TH] 

i3(T^A? 7 ) G F m\ (v \ 2 / ™\ 



«-) (C-7) 



B(T — > ) 4a/27tq; Vx 
where F ~ 1/2 includes QCD corrections and B(T(ls) -> = (2.48 ± 0.06)% [S]. 
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